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Abstract 

We perform the dual transformation of the Yang-Mills theory in d = 3 dimensions 
using the Wilson action on the cubic lattice. The dual lattice is made of tetrahedra 
triangulating a 3-dimensional curved manifold but embedded into a flat 6-dimensional 
space (for the SU{2) gauge group). In the continuum limit the theory can be refor- 
mulated in terms of 6-component gauge-invariant scalar fields having the meaning of 
the external coordinates of the dual lattice sites. These 6-component fields induce a 
metric and a curvature of the 3-dimensional dual colour space. The Yang-Mills theory 
can be identically rewritten as a quantum gravity theory with the Einstein-Hilbert 
action but purely imaginary Newton constant, plus a homogeneous 'matter' term. In- 
terestingly, the theory can be formulated in a gauge-invariant and local form without 
explicit colour degrees of freedom. 



1 Lattice partition function 



Though our objective is the continuum theory we start by formulating the SU(N C ) gauge 
theory on a cubic lattice. The partition function can be written as an integral over all link 
variables being SU(N C ) unitary matrices U with the action being a sum over plaquettes, 

Z(0) = J J] ^imk exp ( £ (3 (Tr C/ plaq + c.c.) / 2Tr 1 j (1) 

links yplaquettes J 

where (3 is the dimensionless inverse coupling. The unitary matrix U p \ aq is a product of four 
link unitary matrices closing a plaquette. 

To get to the continuum limit one writes Uy,^ = exp^iaA^t 01 ) where a is the lattice spacing 
and A a ^t a = is the Yang-Mills gauge potential with t a being the generators of the gauge 
group normalized to Tr t a t b = 5 ab /2, and expands Tr U p i m in the lattice spacing a. As a 
result one gets for a plaquette lying in the (12) plane: 

IV^ + c.o. = \ 
H 2 Tr 1 H \ 2Tr 1 V ') ' V ; 

where F^ v = d^A v — d u A^ — i[A^A v ] is the Yang-Mills field strength. Summing over all 
plaquettes one obtains the partition function of the continuum theory, 



^cont = J DA^exp (-7^2 / d d x Tr 



(3) 



with an obvious relation between the dimensionless lattice coupling (5 and the SU (N c ) gauge 
coupling constant in d dimensions, g\: 

a 4 d 9d 

In this paper we concentrate on the Euclidean £77(2) Yang-Mills theory in d = 3 dimen- 
sions. In this case eq. (f|) reads: 

fi=— 2 - (5) 

The continuum limit of the d = 3 Yang-Mills theory given by the partition function ([[]) is 
obtained as one takes the lattice spacing a — > and (3 — » oo with their product = 4/ (a/?) 
fixed. This quantity provides the theory with a mass scale. It is widely believed (though not 
proven so far) that the theory possesses two fundamental properties: 1) the average of a large 
Wilson loop has an area behaviour with a string tension proportional to gf , 2) correlation 
functions of local operators like decay exponentially at large separations, with a 'mass 
gap' proportional to g\. 

Our aim will be to rewrite the partition function (|1]) in dual variables and to study its 
continuum limit. 



2 



2 Dual transformation 



The general idea is to integrate over link variables Uy,^ in eq. (||) and to make a Fourier 
transformation in the plaquette variables t/ p i aq . This will be made in several steps, one for 
a subsection. 

2.1 Inserting a unity into the partition function 

First of all, one needs to introduce explicitly integration over unitary matrices ascribed to 
the plaquettes, U p \ aq . This is done by inserting a unity for each plaquette into the partition 
function ([]]): 

1 = If / df/ Plaq ^Plaq , W 2 W 4 ) (6) 
plaquettes 

where C/1...4 are the link variables closing into a given plaquette. The 5-function is understood 
with the group-invariant Haar measure. A realization of such a 5-function is given by Wigner 
-D-functions: 

WV)= E (2J + l)D J mim2 (U^)D J m2mi (V). (7) 
J=o,|,i,|,... 

This equation is known as a completeness condition for the D-functions ffl. The main 
properties of the Z)-functions used in this paper are listed in Appendix A. 

Eq. (0) should be understood as follows: if one integrates any function of a unitary matrix 
U with the r.h.s. of eq. (0) over the Haar measure dU one gets the same function but of the 
argument V: 

JdU f(U) 6(U,V) = f(V). (8) 

Using the multiplication law for the Z)-functions (see Appendix A, eq. (|75|) ) one can write 
down the unity to be inserted for each plaquette in the partition function ([I]) as 

1 = / ^Plaq £(2J+ 1) ^ m2 (^la q )^ m3 (^)^4 m4 (^)^ 4m5 (f/3) J D4 mi (^4) (9) 
J J 

where C/1...4 are the corresponding link variables forming the plaquette under consideration. 

2.2 Integrating over plaquette variables 

Integrating over plaquette unitary matrices t/ p i aq becomes now very simple. For each pla- 
quette of the lattice one has factorized integrals of the type 

/ ^Piaq exp ^^P^ + ^^aq j d^u^j = 5 mima jjIi(P)Tj(P), (10) 
where Tj{(3) is the ratio of the modified Bessel functions ||, 
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TjW) 



'2J+1 



W) 



exp 



2J(J+ 1)' 





at /? 



oo. 



(11) 



The quantity Tj(j3) is the 'Fourier transform' of the Wilson action; since in the lattice 
formulation the dynamical variables have the meaning of Euler angles and are therefore 
compact, the Fourier transform depends on discrete values J = 0, 1/2, 1, 3/2, .... However, 
as one approaches the continuum limit (/3 — > oo) the essential values of the plaquette angular 
momenta increase as J ~ ^Jf3 and their discreteness becomes less relevant. Strictly speaking, 
the continuum limit is achieved at plaquette angular momenta J> 1. 

We would like to make a side remark on this occasion. The quantity Tj{(5) gives the 
probability that plaquette momenta J is excited, for given f3. For a typical value used in 
lattice simulations (3 = 2.6 (in 4 dimensions) we find that the probabilities of having plaquette 
excitations with J = 0, 1/2, 1, 3/2 and 2 are 56%, 29%, 11%, 3% and 1%, respectively. It 
means that lattice simulations are actually dealing mainly with J = 0, 1/2 and 1 with a 
tiny admixture of higher excitations. It would be important to understand why and how 
continuum physics is reproduced by lattice simulations despite only such small values of 
plaquette J's are involved. 

We get, thus, for the partition function: 



z= -HP) £ n (2Jp + 1)T 

J Jp plaquettes 

x II / <Wi Dtm 2 {U x ) D% m3 (U 2 ) D% m4 (U 3 ) D J m r 

links I 



- # of plaquettes 



n (2J P + 1) 

plaquettes 



where U 



links I " 

i_4 are link variables forming a plaquette with 



(U 4 

angular momentum J, 



(12) 



'p. 



2.3 Integrating over link variables 

The difficulty in performing integration over link variables in eq. ( |12"D is due to the fact 
that any link enters several plaquettes. In d = 2 dimensions every link is shared by two 
plaquettes, hence one has to calculate integrals of the type 

J dUDi]{U)D J ^) = jj^Sj iJ2 5 kn 5 lm (13) 

for all links on the lattice. We shall consider this case later, in section 4. 

In d = 3 dimensions every link is shared by four plaquettes, hence the integral over link 
variables is of the type 

J dUD^ im2 (U)D^ mA (U)Di% 6 (U)Di- m8 (U) (14) 

where Ja,b,c,d are angular momenta associated with four plaquettes intersecting at a given 
link U, and mi_§ are 'magnetic' quantum numbers, to be contracted inside closed plaquettes. 
In d — 4 dimensions there will be six plaquettes intersecting at a given link but we shall not 
consider this case here. 
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Figure V. 'Even' cubes in checker board order. 



The general strategy in calculating the link integrals (|14j) is (i) to divide by a certain rule 
four D-functions into two pairs and to decompose the pairs of D-functions in terms of single 
D-functions using eq. (j8^ ) of Appendix A, (ii) to integrate the resulting two D-functions 
using eq. (|TjJ) and, finally, (iii) to contract the 'magnetic' indices. Since all 'magnetic' indices 
will be eventually contracted we shall arrive to the partition function written in terms of the 
invariant 3nj symbols. 

There are several different tactics how to divide four D-functions into two pairs, even- 
tually leading to anything from 6j to 18j symbols. In this paper we take a route used in 
refs.p, U, leading to a product of many 6j symbols, although on this route one looses certain 
symmetries, and that causes difficulties later on. The gain, however, is that it is more easy 
to work with 6j symbols than with 12j or 18j symbols. Since important sign factors have 
been omitted in refs. [§, f|] and only final result has been reported there, we feel it necessary 
to give a detailed derivation below. 

In d = 3 dimensions all plaquettes are shared by two adjacent cubes, therefore, it is 
natural to divide all cubes of the lattice into two classes which we shall call 'even' and 'odd', 
and to attribute plaquettes only to even cubes. We shall call the cube even if its left-lower- 
forward corner is a lattice site with even coordinates, (— V) x+V+Z = +1- It will be called 
odd in the opposite case. The even and odd cubes form a 3-dimensional checker board, as 
illustrated in Fig.l, where only even cubes are drawn explicitly. The even cubes touch each 
other through a common edge or link, as do the odd ones among themselves. The even and 
odd cubes have common faces or plaquettes. All plaquettes will be attributed to even cubes 
only: that is the reason for the division of cubes into two classes. 

Let us consider an even cube shown in Fig. 2. 

A, B, C, D, E, F denote its 6 faces, numbers from 1 to 12 denote its links or edges, 
a, b, c, d, e, /, g, h denote its 8 vertices or sites. Correspondingly, we shall denote plaque- 
tte angular momenta by Ja-f, link variables by £/i_i2, and the 'magnetic' numbers of the 
D-functions will carry indices a — h referring to the sites the D-functions are connecting. 

One can write the traces of products of four D-functions over plaquettes in various ways. 
To be systematic we shall adhere to the following rule: Link variables in the plaquette are 
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Figure 2: Elementary even cube 



taken in the anti-clock-wise order, as viewed from the center of the even cube to which the 
given plaquette belongs. If the link goes in the positive direction of the x, y, z axes we ascribe 
the U variable to it; otherwise we ascribe the W variable to it. 

With these rules the six plaquettes of the elementary cube shown in Fig. 2 bring in the 
following six traces of the D-function products: 



Cube 



D 



Jc 

k,- hi-. 



(Ul)D^ k} XU 9 )D^ kg (U 6 )D(X(Ul)} [DfxmDil^DfXiU^D^iUl)} 



D J Z ma (Ul 2 )D J m l md (U 4 )D J m - dmh (Un)D% hm . (Uj)] fe e (U 8 )D J n - ng (U 7 )D J n - nf . 

(15) 

Each link variable £/i_i2 appears in this product twice: once as U, the other time as 
U'. For D(W) we use eq. (|77|) of the Appendix A to write it in terms of D(U). After 
that we can apply the decomposition rule (|82|) of that Appendix to write down pairs of D- 
functions in terms of one D-function and two 3jm symbols. The new D-functions correspond 
to the links and carry angular momenta which we denote by j's. The 3jm symbols have 
'magnetic' indices which get contracted when all indices related to a given corner of the cube 
are assembled together. Though this exercise is straightforward it is rather lengthy, and we 
relegate it to Appendix B. As a result we get the following expression which is identically 
equal to fll5|): 

Cube= ]T (2j 1 + l)...(2 Jl2 + l) 

jl--jl2 

D^ hUg (U 7 )D^_ v M)D%,-^ 



3l J4 3l2 
Oa 'a Z a 

312 js 35 

%e $e 



Jl 39 32 

-o b w b p b 

3a 35 39 

~tf s f w f 



32 3s 3io 

Pc Qc X c 

36 jlO 37 

tg Xg Ug 



34 33 3i\ 

-Td qd Vd 

37 3u 3s 

-u h Vh v h 
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Figure 3: Several cubes combine to produce 6j symbols composed of link momenta j. These 
are the same cubes as in Fig.l. 



x 



J J7 Jll J8 II J6 JlO J7 I I J6 J5 J9 I I Jl2 J8 J5 

[ Je Jf Jd ) \ Jd Jf Jc ) \ Jb Jc Jf J \ Jf Jb Je 

J4 J3 jll 1 f J2 J3 JlO 1 f jl J9 J2 1 f jl J4 Jl2 ) 

Jd Je Ja I I Jd Jc Ja J \ Jc Ja Jb | I </e Jb Ja \ 



(16) 



Here ji-12 are the angular momenta attached to the links of the cube, (...) are 3jm- and 
{...} are 6j-symbols. We see that there is a 6j symbol attached to each corner of the even 
cube; its arguments are three plaquette momenta J and three link momenta j intersecting 
in a given corner. The 3jm symbols involve only link variables j. 

We have, thus, rewritten all twelve pairs of D J -i unctions entering a cube as a product of 
single D J -functions, where j's are the new momenta associated with links. It is understood 
that this procedure should be applied to all even cubes of the lattice. After that, one has 
only two /^-functions of the same link variable U, for all links of the lattice. It becomes, 
therefore, straightforward to integrate over link variables, using eq. (|T3|) . 

It is convenient to integrate simultaneously over six links entering one lattice site, because 
in that way one gets a full contraction over all 'magnetic' numbers. The derivation is, again, 
straightforward but lengthy: the details are given in Appendix C. The result is that the 3jm 
factors in eq. ( p!6|) are contracted with analogous 3jm symbols arising from neighbouring 
even cubes, and produce 6j symbols attached to every lattice site and composed of the six 
link momenta j intersecting at a given lattice site. In notations of Fig. 3 we get for the 
vertices a, b: 





"a" = J1 J4 , "b" = Ji Jy JZ , (17) 



and similarly for other vertices. A sign factor (— l) 2jf should be attributed to every link of the 
lattice. As shown in Appendix C it is actually equivalent to a sign factor (— 1) 2J attributed 
to every lattice plaquette. 
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3 Lattice partition function as a product of 6j symbols 



We summarize here the recipe derived in the previous section. One first divides all 3-cubes 
into two classes, even and odd ones. They form a 3-dimensional checker board depicted in 
Fig.l. All even cubes are characterized by their plaquette momenta J. The edges of even 
cubes have link momenta j; each link is shared by two even cubes. 

To each of the eight corners of an even cube one attributes a Qj symbol of the type 

{ i i i } < i8 > 

where J's are plaquette and j's are link momenta intersecting in a given corner of a cube. 
The rule is that link 1 is perpendicular to plaquette A, link 2 is perpendicular to plaquette 
B and link 3 is perpendicular to plaquette C. Four triades, (JiJbJc), ^JaJo), {.JzJaJb) 
and (J1J2J3) satisfy triangle inequalities. 

To each lattice site one attributes a Qj symbol of the type 



Jl 32 J3 

H 3s je 



(19) 



where j's are the six link momenta entering a given lattice site. The rule is that link 4 is 
a continuation of link 1 lying in the same direction, link 5 is a continuation of link 2 and 
link 6 is a continuation of link 3. Four triades, {3132k), (jijsje), (Jzjda) and (jajds) satisfy 
triangle inequalities. 

Actually, each lattice site has five Qj symbols ascribed to it: four are originating from 
the corners of the even cubes adjacent to the site and are of the type flT8|), and one is of the 
type ©. 

The lattice partition function (|T|) or ( |T2"D can be identically rewritten as a product of the 
Qj symbols described above. Independent summation over all possible plaquette momenta J 
and all possible link momenta j is understood. We write the partition function in a symbolic 
form: 



Z 



I'M 



# of plaquettcs 



E II ( 2 Jp + 1)Tj„ 09) (-1 



^2Jp 



Jp>jl plaquettcs 



n m + 1) 

links 



x n 

even cubes corners 



i; 


3 


3 


} n 1 

' lattice sites 1 


3 


J 


;) 




J 






3 





(20) 



The plaquette weights Tj(j3) are given by eq. ([□]). Apart from the sign factor essentially the 
same expression was given in refs. 0, Q []. The sign factor is equal to ±1 if the total number 
of half-integer plaquettes J's is even (odd). Since plaquettes with half- integer momenta form 
closed surfaces it may seem that the sign factor can be omitted. In a general case, however, 
when one consideres vacuum averages of operators this is not so, therefore, it is preferable 
to keep the sign factor. 



: We are grateful to P.Pobylitsa who has independently derived eq. (EO 
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4 Simple example: d = 2 Yang— Mills 



In a simple exactly soluble case of the 2-dimensional SU(2) theory every link is shared by 
only two plaquettes. Therefore, the link integration is of the type given by eq. fll3"l) : it 
requires that all plaquettes on the lattice have identical momenta J. The partition function 
thus becomes a single sum over the common J: 



Z 



# of plaquettes 

# of plaquettes 



\rpj (/?)]# P la q uettes ^21) 

j 

the number of plaquettes being equal to V/a 2 where V is the full lattice volume (full area 
in this case) and a is the lattice spacing. 

A slightly less trivial exercise is to compute the average of the Wilson loop. Let the 
Wilson loop be in the representation j s . It means that one inserts D^ S (U) for all links 
along the loop. One gets therefore integrals of two D-functions outside and inside the loop, 
and integrals of three D-functions for links along the loop. The first integral says that all 
plaquettes outside the loop are equal to a common J. The second integral says that all 
plaquettes inside the loop are equal to a common J 7 . Integrals along the loop require that 
J, J' and j s satisfy the triangle inequality. We have thus for the average of the Wilson loop 
of area S: 

= E J [T,(/3#E^,[^W/r,W]» 

This is an exact expression for the lattice Wilson loop, however we wish to explore its 
continuum limit. It implies that V/a 2 — > oo, S/a 2 — > oo but S <C V; {3 — > oo, a — > but 
(3a 2 = 4:/g 2 fixed, where g 2 is the physical coupling constant having the dimension of mass 2 , 
see eq. (|J). 

We take V/a 2 — > oo first of all, which requires that only the J = term contributes to 
the sum, with Tq(/3) = 1; consequently all momenta inside the loop are that of the source, 
J' = j s . Taking into account the asymptotics ( |Tl~D of Tj(/3) at large j3 we obtain 



s 



(W 3s (S)} = [T Js ((3)p = exp 



(23) 



which is, of course, the well-known area behaviour of the Wilson loop with the string tension 
proportional to the Casimir eigenvalue. 



5 Dual lattice: tetrahedra and octahedra 

We now turn to the construction of the dual lattice. 



Each 6j symbol of the exact partition function ( 20f) encodes four triangle inequalities 
between the plaquette J's and the link j's. It is therefore natural to represent all 6j symbols 
by tetrahedra whose six edges have lengths equal to the six momenta of a given 6j symbol. 
Four faces of a tetrahedron form four triangles, so that the triangle inequalities for the 
momenta are satisfied automatically. 
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Figure 4: Tetrahedra corresponding to the Qj symbols sitting at vertices a and b. 



h f 




Figure 5: Octahedron dual to the even cube. 

Let us first consider the eight Qj symbols corresponding to the eight corners of an even 
cube. These eight Qj symbols are given explicitly in eq. fljjD with notations shown in Fig. 2. 
Let us represent all of them by tetrahedra of appropriate edge lengths. For example, the 
tetrahedra corresponding to the corners a and b are shown in Fig. 4. We denote the plaquette 
momenta Ja, • •• just by their Latin labels A, B, ... and the link momenta j\,ji, • •• by their 
numerical indices 1,2, ... We notice immediately that the two tetrahedra have a pair of equal 
faces, in this case it is the triangle (A,B,1). Therefore, we can glue the two tetrahedra 
together so that this triangle becomes their common face. The gluing can be done in two 
ways. To be systematic we shall always glue tetrahedra so that their volumes do not overlap. 

In the same way we glue together other tetrahedra. Being glued together the eight 
tetrahedra of the cube form an octahedron shown in Fig. 5. Its center point O is connected 
with six lines to the vertices denoted as A — F; the lengths of these lines are equal to the 
corresponding plaquette momenta Ja-f- The external twelve edges of the octahedron have 
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Figure 6: A tetrahedron corresponding to the lattice site fits precisely into a hole between 
four octahedra corresponding to the four corners of the even cubes adjacent to the site 
(shown in motion). 

lengths equal to the link momenta ji-u- The eight faces of the octahedron correspond to 
the eight vertices of the original even cube. One can say that the octahedron is dual to the 
cube: the faces become vertices and vice versa; the edges remain edges. 

It is clear that in a case of generic J's and j's the octahedron cannot be placed into a 
flat 3-dimensional space. Indeed, we have 6+12 = 18 given momenta, that is fixed lengths, 
but only 7 points defining the octahedron, including the center one. In three dimensions 
that gives 21 d.o.f. from which one has to subtract 3+3 to allow for rigid translations and 
rotations. Therefore, we are left with only 15 d.o.f. instead of the needed 18. [In four 
dimensions the arithmetic would match: 7-4 — 4 — 6 = 18.] 

Each even cube of the original lattice has twelve neighbouring even cubes sharing edges 
with the first one, and with themselves. If we represent the neighbour even cubes by their 
own dual octahedra those will also share common edges. Does this network of octahedra 
cover the space? No, there are holes in between. However, we have not used yet the 6j 
symbols ([19|) made solely of the link momenta j's. If we represent these 6j symbols by 
tetrahedra their triangle faces will coincide with the faces of the octahedra corresponding to 
the even cubes adjacent to the site. For example, if we consider the 6j symbols corresponding 
to the site a (see Fig. 3 and eq. (|T7|)). 



it has a common triangle face (ji, J4, J12) with the octahedron shown in Fig. 5. The other 
faces of this tetrahedron match the octahedra dual to the cubes //, /// and IV, see Fig. 3. 

Octahedra corresponding to the cubes supplemented by tetrahedra corresponding to the 
lattice sites cover the space without holes and therefore serve as a simplicial triangulation, 
see Fig. 6. 

An equivalent view on the dual lattice has been suggested in ref. [[|. One can connect 
centers of neighbour cubes (both even and odd) and ascribe plaquette momenta J's to these 
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Figure 7: Another view on the dual lattice. 



lines. The link momenta j's will be then ascribed to diagonal lines connecting only even 
neighbour sites of that dual lattice, see Fig. 7. 

The dual lattice can be understood in two senses. On one hand, one can build a regular 
cubic dual lattice with additional face diagonals like shown in Figs. 6 and 7, and ascribe J's 
and j's to its edges. On the other hand, since variables living on the links of the dual lattice 
are positive numbers, one can build a lattice with the lengths of edges equal to the appropriate 
angular momenta. We shall always use the dual lattice in this second sense. 

6 Coordinates of the dual lattice as new variables 

In the previous section we have already met with a situation when an octahedron dual to 
a cube did not fit into a 3-dimensional flat space: at least four dimensions were necessary. 
As one enlarges the triangulated complex more dimensions are needed to match the number 
of degrees of freedom. In the limiting case of an infinite lattice one needs 6 flat dimensions. 
This number of dimensions follows from the number of d.o.f. one has to accommodate: at 
each lattice site there are three plaquette momenta J and three link momenta j, and there 
is a one-to-one correspondence between lattice sites and the cubes. 

Therefore, the dual lattice (understood in the second sense, see above) spans a 3- 
dimensional manifold which can be embedded into a 6-dimensional flat space. Notice that it 
is the maximal number of flat dimensions needed to embed a general 3-dimensional rieman- 
nian manifold; it can be counted from the number of components of the metric tensor, which 
is 6 in three dimensions. Only very special configurations of J's and j's would be possible 
to embed into a flat space of less dimensions. 

We are primarily interested in the continuum limit of the lattice theory, that is in the 
small a, large (3 case. It implies that large angular momenta J ~ ^Jfi are involved, and one 
can pass from summation over J's and j's to integration over these variables in the partition 
function (^0|) . We replace 



■7=0,1/2,1,... u 

and similarly for the summation over link momenta j's. 




(24) 
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The next step is to ascribe a 6-dimensional Lorentz scalar field w a (x), a = 1,...,6 to 
the centers of all cubes of the original lattice, see Fig. 7. We shall call them coordinates 
of the dual lattice. They are scalars because in three dimensions the cubes are scalars. 
The argument of the six-component scalar field is the coordinate of the center of the cube 
in question, however, we shall consider w a (x) as continuous functions. Since six functions 
depend only on three coordinates there are three relations between w a (x) at any point; these 
relations define a curved 3-dimensional manifold whose triangulation is given by the set of 
J's and j's. 

We next define six- dimensional angular momenta as differences of w a (x) taken at the 
centers of neighbour cubes: 

2 

J x (x + ^V,^ = w a (x + a,y,z)-w a (x,y,z) = ad x w a + Y d l wa + -- , 

fiz(x + %,y,z+f) = w a (x + a,y,z)-w a (x,y,z + a) = a(d x -d z )w a + 0(a 2 ),(25) 

and so on. The lengths of these 6-vectors are, by construction, the lengths of the edges of 
the dual lattice. 

The six functions w a (x) can be called external coordinates of the manifold; they induce 
a metric tensor of the manifold determined by 

9ij (x) = d iW a d jW a . (26) 
As usual in differential geometry one can define the Christoffel symbol, 

Fi,jk(x) = ^{dj9ik + d k g i:j - diQjk) = diW a djd k w a = (wt ■ w jk ), (27) 
and the Riemann tensor, 

Rijki{x) = -{djd k gu + didigj k — djdigi k — did k gji) + r mjfc r™ — r mj; r^ 

= [(w ik ■ Wji) - g pq (w p ■ w ik ) (w q ■ Wji] -[&<-»■ I]. (28) 
The contravariant tensor is inverse to the covariant one, 

9 ij 9ik = Si (29) 

and can be used to rise indices, and for contractions. The determinant of the metric tensor 
is 

g = det g i:j = ^e ijk e lmn {wi ■ w{) (wj ■ w m ) {w k ■ w n ) , (30) 
and the contravariant metric tensor is 

g ij = l_ e ikl e jmn {wk . ^ . ^ 

There is a useful identity for the antisymmetrized product of two contravariant tensors, valid 
in 3 dimensions: 
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g ik g jl _ g il g jk = ^ ^ / ^ ^ 

The scalar curvature is obtained as a full contraction: 

R = g* k R ljkl = ±(g ik g jl - g il g^ k ) R ljkl 

= —e ijk e i,j ' k ' (w k ■ wy) [2g (w H > ■ w jf ) - e plm e ql ' m ' (w p ■ w iV ){w q ■ w jr ){wi ■ w v ){w m -w m >)\. (33) 
^9 

Recalling that w a is a 6- dimensional vector we can rewrite the scalar curvature in another 
form: 

p ^ , ,ijk ^i'j'k'pll'qmm',,a.,,a',,l3 j3' ')/ Y S .,,8' „,,C„,,C (1A\ 

72^2 /37fe ' ? /3 ' 7 ' <5 ' £r ' W i W i> W j W j' W k W k' W lm W l'm' W p W q- l^ 4 J 

This form makes it clear that the scalar curvature is zero if w a has only three nonzero 
components, which corresponds to a flat 3-dimensional manifold. 

Finally, we would like to point out the Jacobian for the change of integration variables 
from the set of the lengths of the tetrahedra edges, Jf and jf given at all lattice sites, to the 
external coordinates w a . In the continuum limit this Jacobian is quite simple. It is given by 
the determinant of a 6 x 6 matrix composed of the second derivatives: 

J\djf(x) djf(x) = Y\_dw a (x) Jac(w), Jac(w) = det w°j. (35) 

X X 

Since w°j = there are actually six independent second derivatives. The Jacobian is zero 
in the degenerate case when the triangulation by tetrahedra can be embedded in less than 
6 dimensions. 



7 Continuum duality transformation and Bianchi iden- 
tity 

It is instructive at this point to compare the duality transformation on the lattice with that 
in the continuum theory. The continuum partition function (|^) can be written with the help 
of an additional gaussian integration over the 'dual field strength', J?-: 



Z = DJ" DA? exp d 3 x 



93 T 2 



^■i j' 2 i 1 j a (£). A a f) ■ A a -J- c abc A b A c ^ 



(36) 



Eq. (^) is usually called the first-order formalism. 

In the Abelian case when the Ai commutator term is absent, integration over A{ results 
in the ^-function of the Bianchi identity, 



diJij 0, 



or 



(■ijkdiJk 0, 



Jk 



1 



'^ijkJij • 



(37) 



Because of this identity, one can parametrize J k = d k w, and get for the partition function: 
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Z ahc[ = jDw exp J 



d 3 x 



(38 



It represents a theory of a free massless scalar field w. It is in accordance with that in a 3d 
Abelian theory there is only one physical (transverse) polarization. It is easy to check that 
gauge-invariant correlation functions of field strengths coincide with those computed in the 
original formulation. 

In the non-Abelian case integration in A\ is more complicated, and there is no simple 
Bianchi identity for J£ = (l/2)ey However, one can formally perform the Gaussian 
integration over Af H resulting in: 



Z = J DJf det^J- 1 ) exp J d 3 x 



2 v * 



a\2 



- (e tjm d 3 J^) {J-^t ( CHB ^J*) (39) 



where J 1 is the inverse matrix, 



{J' 1 ) t e bcd e klm J d m = 8 ac S u ; detiJ- 1 ) = (det J a k )-\ (40) 

Notice that the second term in the exponent is purely imaginary; the full partition function 
is real because for each configuration J*{x) there exists a configuration with — J^(x), which 
adds a complex conjugate expression. 

We now turn to the discretized version of the dual theory. As explained above, we 
need 6 flat dimensions to embed the dual lattice, and we have introduced 6-dimensional 
momenta J a , see eq. fl2~5]) . These momenta apparently satisfy, e.g., the identity (see Fig.7 
for notations): 

Jz (x, y,z + ^\ - J" (x+ |, y, z^j = w a (x, y,z + a) - w a (x + a, y, z) 

= J? (x + a, y, z + - J* (x + |, y, z + aj , (41) 

and similarly for other components. This is nothing but a discretized version of the Bianchi 
identity, 

e ijk diJ^ = 0, a = l,...,6. (42) 

Therefore, in 6 dimensions one recovers the simple (flat) form of the Bianchi identity for 
the dual field strength. One can say that the complicated (nonlinear) form of the usual 
non-Abelian Bianchi identity is a result of the projection of the flat Bianchi identity onto 
the curved colour space. 



8 Wilson loop 

In this section we present the Wilson loop in the representation j s 



Wj = 



^— Tr P exp i <f d Xi A^T a , (43) 
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in terms of dual variables. 

In terms of the original lattice the Wilson loop corresponds to adding a product of the 
D^ S (U) functions to all links along the loop, with a chain contraction of 'magnetic' indices. 
Because of these insertions, on links containing the loop one has to integrate over three D- 
functions instead of two as for all other links. As a result one gets additional 3jm symbols 
along the loop which combine into the new 9j symbols ascribed to all lattice sites, see 
Appendix D. For example, the 9j symbols ascribed to vertices a and b are (for notations see 
Fig.3): 



"a" 



(44) 



J4 Jl 312 J2 Jl 39 

3i5 Js 3i5 > , "b" = < j' 17 j s jn 

313 3l 3u J I jl8 j'l 3w 

The accompanying sign factors are given in Appendix D. Six triades of the 9j symbol, 
corresponding to all its rows and columns, satisfy the triangle inequalities. 

Contrary to the 6j symbol the 9j symbol cannot be represented by a geometrical figure 
with edges equal to the entries of the 9j symbol. In addition, the link momenta along the 
loop split now into pairs: j\ and j[, ji 5 and j[ 5 , jyj and j[ 7 , and so on. The 'primed' and 
'non-primed' angular momenta satisfy triangle inequalities, with the source j s being the third 
edge of the triangles. If j s is an integer, there is always a contribution with j[ = j\ (and so 
on). If j s is a half- integer one has necessarily j[ ^ j\. 

Thus, there appears to be a fundamental difference between Wilson loops in integer and 
half-integer representations. For integer representations one can proceed as in the vacuum 
case and parametrize the dual lattice sites by the coordinates w a (x) related to angular 



momenta through eq. (|25|). In the half- integer case one cannot uniquely parametrize the 
dual lattice by the coordinates w a (x). In the presence of the Wilson loop in a half- integer 
representation the dual space w a is not simply connected: there is a infinitely thin cylindrical 
'hole' in the dual space along the loop. 



9 Asymptotics of the 6j symbols 

In the continuum limit (3 — *> oo, J, j — > oo one can replace 6j symbols by their asymptotics. 
The asymptotics was ingeniously guessed in a seminal paper by Ponzano and Regge |J and 
later on explicitly derived and improved by Schulten and Gordon 0. The results of these 
works can be summarized as follows. 

First of all one draws a tetrahedron with edges equal to j n + |, where j n are the six 
momenta of a given 6j symbol. It should be stressed that though four momenta triades 
satisfy triangle inequalities, the same triades shifted by | need not. In that case the 6j 
symbol is said to be 'classically forbidden', and it is exponentially suppressed at large j n . 

If j n lie in the 'classically allowed' region, the asymptotics is given by the Ponzano-Regge 
formula: 

I ^ ^ K J y/l2nV(j) In V 2j 4 



(45) 
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Here V(j) is the 3-dimensional volume of the tetrahedron and 9 n is the dihedral angle in the 
tetrahedron, corresponding to the edge j n + ~. Since we are interested in the large- j n limit 
we shall systematically neglect the shifts by |. The tetrahedron volume can be found from 
the Cayley formula: 



288 



5l Si 5l 

51 o j! 51 

51 51 o jf 

51 51 5l o 



(46) 



The dihedral angle corresponding, say, to the edge j\ can be found from 



cos Q\ 



1 5\ + mi -ji-ji-ji- jd + (5i - sm - 5i 
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S(jlj2j3)S(j 1 J 5 ,j (i ) 



(47) 



where 



+ 32 + js) (j2 + js - ji) (ji - 32 + 5s)(ji + 32 ~ 3s) 



(4J 



is the area of the triangle built on the edges ji,2,3- The dihedral angles are defined such that 
< 9 < 7i. Since in section 6 we have defined 6-dimensional angular momenta j a whose 
lengths are the edges of the tetrahedra, we can find the dihedral angles from more simple 
formulae involving scalar products of momenta in the 6-dimensional space. For example, 
eq. (f47D can be rewritten as 



COS6 1 ! 



(ji-5i){5i -5b) - 51(32 ■ 5e) 



\l5hl - (Ji ■ J2) 2 \j3\5l - (ji ■ je) 



(49) 



Notice that the angle is defined to be equal to tx (not 0!) when the two vectors, 3% and je 
coincide; it is zero when they point in the opposite directions. We shall use this formula in 
what follows. 



10 Angle defect 

The Yang-Mills partition function (|20|) is a product of many 6j symbols for each of which 
we use the asymptotic form ( [45|) in approaching the continuum limit. Each cosine can be 
written as a half-sum of exponents of imaginary argument. Therefore, we have to consider 
a sum of a product of many imaginary exponents, 

f[ cos(fi„) = ex P MXX^n) , (50) 

where fl n denotes the argument of the cosine in eq. for the n th 6j symbol, and one has 
to sum over all signs e n — ±1. 

The expression in the exponent of eq. ( |o"0"D can be rearranged as follows: We first pick 
one of the edges of the dual lattice, whose length is a link ji or a plaquette Jp, and combine 



17 



all dihedral angles 9 n related to this edge, as coming from the n th tetrahedron. We then sum 
over all edges of the dual lattice. Therefore, we can write: 

5>„,ft n = ]T J P (j2 e n 9 n (J P )) +J2Ji (E e nOn(ji)) , e n = ±1. (51) 

n P \n=l ) I \n=l ) 

As seen, e.g., from Fig. 7, each plaquette J enters four tetrahedra, therefore the corresponding 
sum over n in eq. (|51|) goes from 1 to 4. Each link j enters six tetrahedra, therefore in this 
case the sum is over six dihedral angles n (j), with appropriate signs e n . 

Let us consider the contribution to eq. ( p0|) when all signs e n = +1, and let us for a 
moment assume that the dual lattice spans a 3-dimensional Euclidean manifold. The sum 
of the dihedral angles about an edge is then equal to Ait — 2n = 2tc in case of summing over 
four tetrahedra, and equal to 6n — 2n = Air in case of summing over six tetrahedra. In the 
first case we get exp(2iiiJ) = (— 1) 2J ; in the second case we get exp(47uj) = (— l) 4j = +1. 
Notice that the sign factor (— 1) 2J compensates exactly the same factor in the partition 
function fl2Cf ). We conclude that, if the configuration of the momenta is 'flat', there exists a 
contribution to the sum fl5"0|) that does not oscillate with varying J's and j's. In fact, there 
are exactly two such contributions corresponding to taking all signs e„ = ±1 simultaneously. 
Contributions of any other choice of the signs are oscillating fast at large J's and j's, and 
thus die out in the continuum limit. 

A generic configuration of momenta cannot be embedded into a flat 3-dimensional space, 
however. Therefore, the sum of dihedral angles about the edges J and j will, generally, differ 
from 27r and 4n, respectively. These differences are sometimes called angle deficiencies or 
angle defects (we shall use the second term). Let us denote them: 



®(J) = E^( J )- 27r > ( 52 ) 

n=l 

©0') = E0nO')-47T. (53) 



n=l 



Our task is to point out contributions to eq. (|50D that survive the continuum limit in a 
general case when the dual lattice is a curved 3-dimensional manifold. To be more precise, 
we have to consider the sum of all momenta on the lattice times their angle defects, 



exp % 



(54) 



and to find the contribution of the order of a 3 to this exponent, where a is the lattice 
spacing. The 0(a 3 ) order is needed to compensate for the 1/a 3 factor arising as one goes 
from summation over the lattice points to integration over the 3-dimensional space. 

In the continuum limit we assume that the momenta are given by the gradients of a 6- 
component function w a (x) having the meaning of the 6-dimensional coordinates of the dual 
lattice sites, see eq. (^). If we restrict ourselves to the first terms in the gradient expansion 
in eq. (p5|) , the momenta will be expressed only through three vectors, d x w a , d y w a and 
d z w a . Three vectors define a flat 3-dimensional space; therefore, the angle defects G are 
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zero in the first-derivative approximation. To get a non-zero angle defect it is necessary to 
expand the momenta in eq. (p5|) up to the second derivatives of w a . We shall see that it is 
also sufficient in three dimensions. 

Since the angle defects 0's vanish if j's are taken to the first approximation of the gradient 
expansion, it means that the expansion of G's starts from terms linear in the lattice spacing 
a. According to eq. (25) the expansion of the momenta also starts from terms linear in 
a. Therefore, one can expect that the expansion of the exponent in eq. (0) starts from 
the 0(a 2 ) terms. Were that so, the configuration would be too 'ultraviolet' and would not 
survive the continuum limit. Fortunately, there appears to be an exact cancellation of all 
0(a 2 ) terms in the sum over several neighbour edges of the dual lattice, so that the exponent 
in eq. ( |54j ) proves to be finite in the continuum limit. 

We next embark a rather tedious enterprise of calculating the angle defects about six 
plaquette J's in a cube (each entering four tetrahedra), and about twelve link j's being 
edges of that cube (each involved in six tetrahedra, see section 5). Unluckily, it seems that 
it is the minimal elementary group that is being repeated through the lattice. It means 
that we have to compute as much as 6 • 4 + 12 • 6 = 96 dihedral angles, expressing them 
through the first and second gradients of the 6-component function w a using eqs. ( P5j , W). 
This formidable calculation has been performed by heavily exploiting Mathematica. The 
intermediate results are very lengthy and we do not present them here. However, the final 
result is beautiful. From a direct calculation we obtain: 



exp i 



p i 



exp i 



points x 



exp - J d 3 x yjg{w) R{w), 



(55) 



where g is the determinant of the induced metric tensor as given by eq. fl5U|), and R is the 
corresponding scalar curvature given by eq. (|33|). Actually, we obtain the expression for the 
l.h.s. of eq. (|55f) in the form of eq. (|3~3"D (written in components, 384 terms!) from where we 
recognize that we are dealing with the scalar curvature. 

In fact this result is a concrete realization of a more general theory developed many years 
ago by Regge 0, In these papers it was shown that the l.h.s. of eq. fl55| ) should be equal 
to its r.h.s. for any simplicial triangulation, provided it has a smooth continuum limit. No 
relation of the scalar curvature R to any concrete triangulation was given, though. We feel 
that it is the first time that this ingenious relation has been derived explicitly for a concrete 
triangulation, and the continuum limit shown to exist. 



11 Full partition function 



Having dealt with the 6j symbols of the partition function ( P0|) we now turn to the weight 
factors Tj((3). According to eq. (|TT| ) at large 8 and J we have: 



II Tj(J3) = exp 

plaquettes 



V — 

^ 3 

plaquettes " 
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exp 



d 3 x 2 J[ 

3 







exp 



2 

Jd 3 x^diW a diW a + 0{a 2 



(56) 



where the relation (|5|) between (3 and the physical coupling constant g\ has been used, 



together with the gradient expansion for the angular momenta (|25|). Combining eqs. 
56) and using diW a diW a = gu(w) we get finally for the Yang-Mills partition function: 



Z = J Dw a (x) Jac(u>) g(w) 4 exp J d 3 x 



(57) 



The second term is the Einstein-Hilbert action with a purely imaginary Newton constant; 
it is invariant under global 6-dimensional rotations of the external coordinates w a (x) and, 
more important, under local 3-dimensional diffeomorphisms w a (x) — > w a (x'(x)). 

The first term in eq. ( |5TD can be viewed as a 'matter' source, 



gi 



d 3 x gi 



gi 



d 3 xJg-T i ig 



i j ■ 



(5? 



2 J 2 

with the stress-energy tensor T^y/g = 5^ violating the invariance under diffeomorphisms. 
Since it is homogeneous in space it can be called the 'ether'. 

The functional measure in eq. ([5j]) arises from two sources. One factor is the Jacobian 
for the change of variables from the tetrahedra edges J's and j's to w a , see eq. fl3~5|). The 
other factor arises from the tetrahedra volumes in the asymptotics of the 6j symbols (flBl). 
In the continuum limit the tetrahedron volume can be written as V(j) ~ ^g, and there are 
5 tetrahedra per lattice site, see section 3. 

Once the partition function is written in covariant terms one can forget the origin of 
the external coordinates w a (as the coordinates of the dual lattice) and consider the metric 
tensor as independent dynamical variables over which one integrates in eq. (0). The 
Jacobian for this change of variables can be easily worked out: in fact it is the inverse of 
Jac(u>) introduced in eq. (p5|). As a result we get the integration measure for the partition 
function (| 



Dgu g 



instead of 



Dgu g~ 



(59) 



which would be the invariant measure in 3d. We shall get an independent check of the power 
— | in the next section. However, it is anyhow a local counterterm not affecting the physics. 

We stress that the partition function written in terms of the metric tensor does not 
contain explicit colour degrees of freedom. Nevertheless, implicitly the theory does 
contain three gluons at short distances. 

Indeed, let us make a simple dimensional analysis of eq. (|57|). The dimension of the first 
term in eq. (|57D is g\d 2 w 2 (we are just counting the number of derivatives and the overall 
power of w)\ the dimension of the second term is d 3 w l . At short distances where quantum 
fluctuations of w a (x) vary fast, the second term dominates the first one. Meanwhile, the 
second term is a fast-oscillating functional at nonzero R. Therefore, the leading contribution 
to the functional integral arises from zero-curvature fluctuations of w a , that is essentially 
from the 3-dimensional w a . Being plugged into the first term, the three components of 
w a describe three massless scalar fields. These fields correspond to three gluons of SU(2) 
with one physical (transverse) polarization. It should be paralleled to eq. (|38D for free 
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electrodynamics. This is the correct result for the non-Abelian theory at short distances in 
three dimensions. 

At large distances or at low field momenta the dominant term is, on the contrary, the 
first one as it has less derivatives. It describes six (instead of three) massless scalar degrees 
of freedom. It is the correct number of gauge-invariant degrees of freedom in the SU(2) 
theory. However, the theory remains strongly nonlinear, and it is not clear so far whether 
massless modes survive in the physical spectrum. 

12 Quantum gravity from first-order continuum for- 
malism 

In this section we give another derivation of the partition function fl57|) directly in the 
continuum theory starting from the first-order formalism, see section 7. We shall show that 
the two terms in the exponent of eq. (|36|) are in fact in one-to-one correspondence with the 
two terms in eq. fl5T|), and that the integration measure coincides with that of eq. (|59D. 

Actually, it has been already derived in the previous section that the first terms of 
eqs. (|36|) and (|57|) are equal: 

2 2 2 

S 1 = -^jd 3 x ( J? f = -jJ d 3 x d lW a d lW a = ~jj d 3 x 9ii . (60) 
Let us derive a less trivial relation for the second terms: 

S 2 = l -J d 3 x e^ k Jt(d 3 A a k - d k A) + e abc A b A c k ) = \\d z x y/gH (61) 

This derivation will be done in two steps. We shall first show, following Witten JJ, that 
the l.h.s. of eq. (|6T|) can be presented as a certain Chern-Simons term. Second, we shall 
show that it is formally equal to the Einstein-Hilbert action. A subtle question about the 
integration measure will be discussed at the end of the section. 

The l.h.s. of eq. ( |BTD is apparently invariant under ordinary gauge transformations: 

8 At = -d^u" + e abc J>A\ = -Df(A)u b , 5 J? = e abc u b J^ (62) 

where Df b (A) = diS ab + t acb A\ is the covariant derivative. 

Less evident, it is also invariant under the following local transformation: 

5 Jt = -d lP a - e abc p b A c l , 5A$ = 0. (63) 
Indeed, after integrating by parts we obtain the following variation of the action: 

6S 2 = ~ J d 3 x p a e l]k Df(A)F b k , F b k = djA* - d k A b + e bcd A)A d k . (64) 

This variation is zero owing to the Bianchi identity, e^ k Df b F b k = 0. 

The two transformations combined form a 6-parameter gauged Poincare group, called 
I SO (3). Indeed, let us introduce three 'momenta' generators Pi and three 'angular momenta' 
generators Li satisfying the Poincare algebra, 
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[P a Pb] = 0, [L a L b ] = ie abc L c , [L a P b ] = ie abc P c . (65) 
We next introduce a 6-component vector field Bf 

B, = J^ + A°L^B?T°, Z-={£ a T 3 l = atl! A «*) 
Its gauge transformation has the standard form: 

Bi -> S^B.S + zS^S, 5 = exp[tp a P a + iu a L a \. (67) 
Using the Poincare algebra (EB) it is easy to check that its infinitesimal form coincides with 



eqs. flU Eg) 



Since the l.h.s. of eq. (pl|) is invariant under these 6-parameter transformations, it can be 
rewritten in an explicitly /SO (3)-invariant form. To that end we notice that the invariant 
tensor of this group is 




where "1" is a unit 3x3 matrix. This matrix defines a scalar product, B a M a ^C f3 , which is 
invariant under global (x-independent) transformations (|67|). With the help of this invariant 



tensor we build a local gauge-invariant action having the form of the Chern-Simons term: 

S 2 = l -J d\ e^ k M aP B? (dX + \fP s B]b£} , (69) 
where Fp = —Ffy are the /SO (3) structure constants. Explicitly, 

pa n TP a c 

r bc ~ u 5 r 3+b,c ~ fc abc 

^3+b,c = -^3+6,3+c = 0) -^3+6,3+c = e abc- (70) 



Using the definition fl6"6]) it is easy to check that eq. (|69| ) coincides with the l.h.s. of eq. (pip, 
however it is explicitly invariant under the 6-parameter gauge transformation (|67|) . 

Eq. ( |B"9"D has the form of the Chern-Simons term in a Yang-Mills theory. Though our 
derivation above is for the gauge group SU(2) it is trivially generalized to any Lie group: 
to that end it is sufficient to replace the SU(2) structure constants e a bc by the structure 
constants f a b c of the gauge group under consideration. We also note in passing that in 
four dimensions the mixed iJ^ u F^ u (A) term of the first-order formalism also posseses an 
additional local symmetry. To unveil it, it is sufficient to replace the scalar parameter p a 
in the transformation ( |63"D by a 4- vector parameter p a : the invariance is again due to the 
Bianchi identity, this time in four dimensions. 

The second step in the derivation is more standard. Introducing the dreibein e", a = 
1, 2, 3, satisfying the condition e^e bl = 5 ab , so that the metric tensor is = e?e", and the 
connection 

< b = \z a \^\ ~ d k e\) - le bk (d ie a k - d k e1) - \e ak e bl - d ie c k ), (71) 
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one can identically rewrite y/gR as 



JgR = el {drf - d k u] + e abc Jfr%) ( 72 ) 

where u a i = w° = \e a bcOj\ c . Finally, one notices that, if one makes an identification of the 
dreibein with the dual field strength, e" = Jf, and of the connection with the Yang-Mills 



potential, uif = A 1 }, then eq. ([72]) takes exactly the form of the l.h.s. of eq. (|61|). This parallel 
has been first noticed in ref. ||10|| . 

There is a subtle point in this formal derivation, however. The use of the first-order 
formalism implies that one integrates both over Jf and over A\ (see eq. (|36|) ) or, equivalently, 
over the dreibein and over the connection, independently. Meanwhile, the use of the Einstein 
quantum gravity implies that the connection is rigidly related to the dreibein via eq. (f7lD , 
moreover, we have explicitly used this relation in the above derivation. In ref. || Witten has 
presented arguments that one can, nevertheless, integrate over the connection as independent 
variable. However, the arguments rely upon the use of the equations of motion (one of which 
is the relation (|TTD), and that might be dangerous in full quantum field theory. 

The present paper gives a different kind of argument that the two approaches are in fact 
equivalent. We start with the Yang-Mills partition function. On one hand it can be presented 
in the first-order formalism where one integrates independently over Jf (the dreibein) and 
over Af (the connection). On the other hand we have shown that the Yang-Mills theory 
is equivalent to quantum gravity where one integrates over the external coordinates w a , or 
over the metric tensor, or over the dreibein only. 

Since pure gravity can be rewritten as a Chern-Simons term (|69|), it is actually a topo- 
logical field theory ||, with no real propagating particles. It is the 'ether' term that violates 
the invariance under diffeomorphisms and restores the propagation of gluons, as it should 
be in the Yang-Mills theory, see the end of the previous section. 

Finally, we would like to remark that the integration measure ([59]) could be anticipated 
from the first-order formalism as well. Indeed, integrating in eq. (|39| ) over A® one gets 
eq. (]iO"D, where the integration measure over the dreibein is (det J)~ 3 ^ 2 ~ g~ 3 ^- The Jaco- 

ij ■ 



bian for the change of variables from the dreibein to the metric tensor is del ~ &9ij 9 ^ 2 



Adding the powers we obtain: — | — ~ = — ~, as in eq. 



13 Conclusions and outlook 

In this paper, we have studied the dual transformation of the SU(2) Yang-Mills theory in 3 
dimensions, both from the continuum and lattice points of view. 

On the lattice, one can introduce dual variables being the angular momenta of the pla- 
quettes (J's) supplemented by those associated with the links (j's). The partition function 
can be identically rewritten as a product of 6j symbols made of those angular momenta. A 
Wilson loop corresponds to taking a product of 9j symbols replacing the 6j symbols along 
the loop. One can construct a dual lattice made of tetrahedra whose edges have the lengths 
equal to J's and j's; the tetrahedra span a 3d curved manifold which can be embedded into 
a flat 6d space. 

In the continuum limit the angular momenta are large, and we have introduced continuum 
6d Euclidean external coordinates w a (x) to describe the curved dual space. The Bianchi 
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condition for the Yang-Mills field strength has been shown to be trivially soluble in flat six 
dimensions. 

At large angular momenta one can use the asymptotics of the 6j symbols, given by 
Ponzano and Regge. Using a specific simplicial triangulation of the dual space (as dictated 
by the original lattice) we have shown that the product of the 6j symbols does have a 
smooth continuum limit which appears to be the Einstein-Hilbert action, with the metric 
tensor and the scalar curvature R expressed through the flat external coordinates w a (x). 
This result cannot be considered as particularly new (it is the cornerstone of the Regge's 
simplicial gravity), however, to our best knowledge it is the first time that the result has 
been explicitly derived from a concrete triangulation of the curved space, and the continuum 
limit shown to exist. We have also found the integration measure for the continuum limit. 

The continuum Yang-Mills partition function can be rewritten as a quantum gravity 
theory but with an 'ether' term violating the invariance in respect to general coordinate 
transformations or diffeomorphisms. This term, however, revives gluons at short distances, 
in contrast to the topological pure gravity theory where no particles propagate. 

The presentation of the Yang-Mills theory in a quantum gravity form ([57]) is explicitly 
colour gauge-invariant since the metric tensor of the dual space is colour-neutral. We have, 
thus, formulated the Yang-Mills theory solely in terms of colourless 'glueball' degrees of 
freedom |. It turns out to be an interacting theory of six massless scalar fields. Nevertheless, 
at small distances it correctly reproduces the propagation of gluons. It is not clear to us 
at the moment how to proceed best in order to reveal its large- distance behaviour. Let us 
indicate a few possibilities. 

One possibility is to exploit the fact that the pure quantum gravity theory is topological, 
therefore essentially a free theory. One can try to make a perturbative expansion in g% about 
it. 

Another possibility is to make use of the fact that the Chern-Simons term ( |55D can 
be obtained from integrating over heavy fermions, in this case belonging to some IS 0(3) 
representation. The subsequent integration over bosonic fields Ai, Ji is trivial since there is 
no kinetic energy term for those fields: the result would be a local four-fermion theory with 
infinitely heavy fermions; it might be soluble, at least in the large- N c limit. 

Probably the most promising possibility is to pursue the analogy with and methods of 
quantum gravity. One can average eq. (57) over 3d diffeomorphisms: the second term is 
invariant, the first term is not. Integrating the first term over diffeomorphisms produces 
diffeomorphism-invariant effective action containing growing powers in the curvature. The 
effective action may lead to a nonzero v.e.v. of the scalar curvature, and that may yield 
a mass gap for the diffeomorphism-non-invariant correlation functions, like the correlation 
functions of F^ u . 

There are several other tasks for the future, lying on the surface. First, it would be 
interesting to generalize the present approach to colour groups other than SU(2). In view of 
the sad fact that the theory of the "6j symbols" for higher Lie groups is not too developed it 
will be probably difficult to make a straightforward generalization of the lattice formulation. 
A more promising approach would be to start from the first-order formalism, the more so 
that the wide local symmetry revealed in section 12 can be directly generalized to any Lie 



2 A somewhat similar line was developed in ref. [[IT] for the 3+1 dimensional Yang-Mills theory in the 
Hamiltonian approach; see also ref. [O] . 
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group. Second, it would be interesting to make a transformation similar to that of this paper 
in d = 4. The lattice 6j symbols have been known for a while in this case [|J (for the SU(2) 
colour), however it again seems that the first-order formalism is a more promising start, due 
to the additional gauge symmetry noticed in section 12. 

We are grateful to Pavel Pobylitsa for many fruitful discussions. D.D. acknowledges a very 
useful conversation with Ben Mottelson. V.P. is grateful to NORDITA for the hospitality 
extended to him in Copenhagen, in particular for a partial support by a Nordic Project 
grant. The work was supported in part by the Russian Foundation for Basic Research grant 
97-27-15L. 



Appendix A. D-functions, 3jm, 6j and 9j symbols 

Wigner D-functions are eigenfunctions of the square of the angular momentum operator 
(written in terms of, say, three Euler angles a, (3, 7), 

J 2 D J mn (a,(3^) = J{J+\)DiJa,P,i), J = 0, ^ 1, |., -J< m , n <+J, (73) 

and can be said to be eigenfunctions of a spherical top; they are (2 J + l) 2 -fold degenerate. 
The 'magnetic' quantum numbers m, n have the meaning of the projections of the angular 
momentum of a spherical top on the third axes in the 'body-fixed' and 'lab' frames. One 
can parametrize a 2 x 2 unitary matrix by Euler angles as 

U = exp(zar 3 ) exp(i/3r 2 ) exp(«7r 3 ). (74) 

It is convenient to use the unitary matrix U as a formal argument of the D-functions. Their 
main properties are: 

• Multiplication law: 

D^UxUz) = D'l m {Ui)D^ d {U2) (summation over repeated indices understood). 

(75) 

• Unitarity: 

D J M (U ] ) = (p J lk {U)) * ( " * " denotes complex conjugate). (76) 

• Phase condition: 

(DfkiU))* = (-l) l - k D\_ k (U), D J kl (l) = 5% J+1 \ (77) 

• Orthogonality and normalization: 

J dUD(l(tf)D% n (U) = 2^n:^A **n Sim- (78) 
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Integration here is over the Haar measure: 

J dU... = J d{SU)... = J d(US)...; JdU = l. (79) 

• Completeness (the 5-function is understood in the Haar measure sense): 

m V) = £(2J + l)D J kl {tf)Di{V). (80) 
j 

• Matrix element: 

JdUD* l( U )D^ ( U )D^(U) = ( * * * ) ( £ ll * ) , (81) 

where (...) denote 3jm symbols. 

• Decomposition of a direct product of ireps: 

(82) 

The last two factors may be replaced by D J _ d _ c {U^) using eq. (|77|) . 

The 3jm symbols are symmetric under cyclic permutations of the columns. An inter- 
change of two columns gives a sign factor: 

( J : J ; ■>'• ) = (-1)*+*+* ( J ; j ; J: '' ) , etc. (83) 
\ k I m I \ I km) 

If one changes the signs of all 'magnetic' quantum numbers or projections, the 3jm symbol 
also gets a sign factor: 

( t ? t ) - ( \ -i t ) ■ < 84 > 

A "practical" definition of the 6j symbol {...} is via a contraction over projections in 
three 3jm symbols: 



E(-i) 



kirn 



i4 _ fc+i5 _; +i6 _ m J 5 jl J6 W J6 32 34 \ I 3l 33 3b 

I / p — m J I m q —k J I k r —I 

= ( 3\ 32 33 \ f ji 32 33 j 

V ~P -Q ~r ) 1 J4 js 36 J ' 



(85) 



The summation over projections k, I, m is such that p = m — I, q = k — m and r = I — k are 
kept fixed. 

Another definition of the 6j symbol is via the full contraction of projections in four 3jm 
symbols: 



26 



^_iy4+n+j5+o+j 6 +p 

klmnop 





(86) 

Since the three j's of any 3jm symbol satisfy the triangle inequalities, e.g. \j\ — J2I < 
jz < h +32) etc, the following four triades of the 6j symbols have to satisfy the triangle 
inequalities: (jijtfa), (Jdde), (Jdda) and (Jajds)] otherwise, the 6j symbol is zero. 

The 6j symbols are symmetric under permutation of any of two columns and under 
interchange of the upper and lower arguments simultaneously in any two columns, e.g., 



I Jl 32 33 I _ I Jl 33 32 I _ I J4 32 36 

\ 34 h 3e J \ 34 3e 35 J \ 3i 35 h 
A full contraction of six 3jm symbols yields the 9j symbol: 



etc. (87) 




y ( Ji 32 33 \ I 34 Js 36 \ I 37 3s 39 \ I 3i 34 3i 
^ \ k I m J \ n o p J \ q r s I \ k n q 

jl 32 33 
34 35 36 
37 38 39 

9j symbol is symmetric under transposition and under even permutations of rows and 
columns; under odd permutations it acquires a sign factor (— l)i 1+ ---+J9. As follows from 
the definition, six momenta triades corresponding to the rows and columns of the 9j symbol 
satisfy triangle inequalities. 

A convenient reference book on D-functions, 3 jm, 6j and 9j symbols is ref . [1J from where 
we have borrowed the definitions. 



Appendix B. 6j symbols in an 'even' cube 

In this Appendix we make the decomposition of two plaquettes D J -functions into a sum of 
single D J -functions labelled by link angular momenta j. Then we assemble the arising 3jm 
symbols into 6j symbols attached to the corners of the even cubes. The notations are given 
in Fig.2. 

We find it convenient (though not necessary) to write the decomposition for the pairs 
containing £71,4,12,6,7,10 (these are links sitting at lower left and upper right corners of the 
cube) in terms of D(U), and the rest in terms of D{U^). 

Exploiting eq. ( [82] ) of Appendix A we get: 



-iy--* £(2j! + 1) ( _^ J . A J _ B . )( _^ J . A J _i )(-ir-^Dii 0b (u 1J 
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=(-D^i;(% + i)(4 j : 4) (4 j c 4H-^ 

=<-i)~p* + i>(4 t 4) (4 t D 4K-<^ 

^-„,(c/i) Ai(£4) 

= (-!)»-», + 1) ( 4 J * ) ( 4 J _* ) ^,-.,(^), 

(-ir^B%+i)(4 £ 4)(4 t 4 ) (-1)"-"^™. 



J6 

J7 ■Tf Jd \ ( J7 Jf Jl 

.17 



-(-^- l ^^)[Z h Z X){-u 9 X -?J(- 1 )"^u(^ 

D J Z nh {Us) D J m l me (Ul) 

= (-!)«- V(2j 8 + 1)( Js Jf Je )( j8 Jf Je W* 04), 
y ! JS J \ -v e n e -m e J \ -v h n h -m h J Vh ' vA 8h 

=(-^n^){t b t -!)(-; % -° f )^M), 

= (-l)^E(2j'io + l)f h l \ D J l){ jW x ^ Jc k )(-^DZ 9 (U W ), 

j 10 \ X c l c ft c / \ Xg lg Kg j 

D J m E dmh (Un) Df h l(Ul) 

= ( _l)^V(2j 11 + l)( ju Je J ° )( 311 Je J f W" -niUu), 
K ' V ~Vd m d ~h j V ~Vh m h -k j Vh ' ydK ' 
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-i)— £(2^+1) ( J _i f a ± a ) ( J _i J » X ) (-ir^D-uuu) 



.112 



(89) 

We now combine together 3jm symbols related to the same vertices (they are marked by 
appropriate indices of the projections a, b, c, d, r, f, g, h), three 3jm symbols for each vertex, 
together with appropriate sign factors. The three 3jm symbols per vertex combine into Qj 
symbols, one for each vertex of the cube. 

Vertex a 



Related to vertex a are the factors 

^_^ya+ja+rn a -o a -r a -Za 



ji J a Jb \ ( U j e J a \ I 3i2 Jb Je 
-o a i a ~ja I \ ~T a m a -i a I { -z a j a -m 



[we use o a + r a + z a = 0, make cyclic permutations in all 3jm symbols, and change the 
summation indices i, j, m — > —i, —j, —m] 



ma ( Jb ji Ja \ f Ja H Je \ ( Je ju Jb \ 

V Ja -O a ~ia ) \ «'a ~T a ~™a j \ ™ a ~Z a -j a J 

.1\-Ja-Jb-Je ( h k 312 j j jl J4 jl2 1 _ /gg\ 

\ ®a ' a %a I 1 Je Jb Ja J 



In the last transformation the definition of the Qj symbol (j85|) has been used. 
Vertex b 

Related to vertex b are the factors 



ib,3b,kb 



( h Ja Jb \ ( 32 Ja Jc \( jg Jc Jb 
o b =ib-j b \ -o b % h -jb ) I ~Pb ib -h / 1 ~Wb h -jb 



[we interchange the first two columns in the first 3jm symbol and change the signs of all its 
projections; it doesn't change the sign of the 3jm's. Also, we make cyclic permutations of 
the last two 3jm symbols, and change the summation indices i, j, k — > —i, — j, —k] 

= J2(-l)~ i+j ~ k [insert 1 = (-l) 2 - 7 ^'" 

( Ja ji Jb \ ( Jb j\ Jc \ ( Jc jg Ja \ 
\ i o b -3 ) I 3 ~Wb -k ) \ k -p b -i ) 



)Js _j a _j c I 31 39 32 \ I 3i 39 32 
I -o b w b Pb \ Jc Ja Jb 



(91) 



In each case we combine the three 3jm symbols and the sign factors so that they suit 
the definition of the 6j symbol given in Appendix A, eq. (^). 

An important property of the sign factors is the following: if ji, Ja, Jb enter one 3jm 
symbol, there is an equality: 
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(_ 1 )±2j 1 ±2J A ±2J B = 1; ^ 

where all signs are possible. This is because out of three momenta either zero or two moments 
are half-integer. Another important property is that, if J is the momentum entering a certain 
3jm symbol, and m is its projection, then (— i) 2J±2m = This is because J and m are 
either integer or half-integer, but simultaneously. 

Below we cite without detailed derivation (which is quite similar to those above) the 
expressions for other vertices of the cube. 

Vertex c 



Vertex d 



Vertex e 



Vertex / 



Vertex g 



-1) 



Ja+Jd-Jc 



J a—Jd—Je 



Je—Jb—Jf 



-1) 



Jb+Jc—Jf 



-1 



.Jc+Jd+Jf 



J2 J3 3l0 
Pc Qc -^c 



J4 33 Jll 

-r d Qd Vd 



312 38 35 

%e V e S e 



3d 35 39 
-tf Sf w f 



36 JlO 37 

tg Xg Ug 



32 33 3w 
Jd Jc Ja 



J4 33 3n 

Jd Je Ja 



312 38 35 

Jf Jb Je 



36 35 39 

Jb Jc Jf 



3e J 10 37 
Jd Jf Jc 



(93) 



(94) 



(95) 



(96) 



(97) 



Vertex h 



- (_\\Je+Jf-Jd 37 311 38 I 37 311 38 I /g^ 
" "' \ ~Uh Vh v h J \ J E J F J D J ' 

Combining all these factors we get eq. ( |I6"D corresponding to the cube. 
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Appendix C. 6j symbols at the lattice sites 



In this Appendix we show how integration over link variables in eq. (|T(J) combine, together 
with the 3jm factors, into 6j symbols composed of the link momenta j, one for each site of 
the lattice. The notations are given in Fig. 3. 

Let us consider integration over link variables ^71,4,12,13,14,15 entering the vertex a shown 
in Fig. 3. This vertex is an intersection of four even cubes denoted in Fig. 3 as I, II, III and 
IV. Link 1 is common to the cubes I and II, link 4 is common to I and IV, and so on. 

The analytical expression for the cube I is given by eq. (|i~6|). The factors relevant to 
vertex a are 




(99) 



It is not necessary to compute anew corresponding expressions for the cubes II — IV. 
It is sufficient to draw a correspondence between the links and the sites of other cubes with 
those of the cube I. For example, link 1, as seen from the viewpoint of cube II, is analogous 
to link 7 of cube 7; the vertex a from the viewpoint of cube 77 is analogous to vertex h of 
cube 7, and vertex b is analogous to vertex g. In the table below we give the list of the 
'analogs' of links in cubes 77 — IV to those of the cube 7. 



II 


I 


1 


7 


13 


11 


14 


8 


a 


h 



III 


I 


12 


10 


14 


2 


15 


3 


a 


c 



IV 


I 


4 


6 


13 


9 


15 


5 


a 


f 



Having this table of correspondence we can immediately read off from eq. 
pressions relevant to the vertex a, arising from the cubes 77 — IV: 



the ex- 



from cube 77 : 



,313 

Va 



-vp 



-u a 
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from cube 777 : ^({/^^.^(f/i^n^i^) 



3l2 J 14 Jl5 
%a Pa Qa 



(101) 



from cube IV: D^ Ub (U 4 )D^_ S£ (U^D^_ w ^ ( ^ f ^ 



(102) 



Integrating over [7^4. 



12,13,14,15 



we get: 



(103) 
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/ dU^Dt d (U,)Dix(U 4 ) = ^-(-i^Sr^Sr,,-^ (104) 

/ dU 12 DZ e (U 12 )D^ Xe (U 12 ) = ^^(-iy*- x «6 Za ,- Xa 5 Ze ,- Xe , (105) 

6 

dUnD^ yr M^ D -^- W Mz) = ^f^i-^-^^-yJ^-y,, (106) 

I dU u D^_ Va (Ul)D^_ Pa (Ul) = -J2^(-iy^ 6pa> _ v J pa> _ Va , (107) 

The four 3jm symbols in eqs.(^-^U^) get now fully contracted over all indices. This 
results in a 6j symbol according to eq. ([86]) of Appendix A. Indeed we have for vertex a: 



_ ^ ^ o-\-t-\-z q v y 

orqvyz 



Jl J4 Jl2 \ / Jl Jl3 Jl4 \ I 312 3U 315 \ I J4 Jl5 Jl3 

or z / I o y v J \ z —v q I 1 r — g — y 



(109) 



[we note that o + r + z = 0; we change the summation variable y — > — y, and interchange the 
last two columns in the second 3jm symbol and the first two columns in the last two 3jm 
symbols - that gives sign factors ^ (_i)iia+ii4+ius an d (_i)i4+ii3+ii 5 _ Finally, 

we insert two unities in the form of 1 = (— l) 2v ~ 2 ^ li and 1 = (— 1) 2 9 _2 J' 15 ] 

— (_]^^jl+j4+jl2+il3-il4-jl5 ^ 

( jl H 312 \ ( jl jlA jl3 \ ( jib jl jl3 \ ( jl5 jlA jl2 \ 

I o r z J y o v —y J y —q r y J \ q —v z J 

= (_^j'i+i4+ii2+ji3-ii4-ii5 J 3i J4 3i2 I (110) 

\ Jl5 Jl4 jl3 J 

[since ji 2 , ju and j'15 came from one 3jm symbol one can use the equation (see eq. fl92j)) 

(_l)il2-il4-il5 _ (_^-il2+jl4+jl5j 

= (_\\h+h-3i2+hz+h4+hs J 3i 3i J12 I (111) 

\ Jl5 Jl4 jl3 J 

This is the final result for the vertex a: the six angular momenta ascribed to the six links 
entering this vertex combine to produce a 6j symbol. 

Similarly, one can treat the vertex b, see Fig. 3. Links labelled by numbers 1, 2, 9, 16, 17, 18 
enter this vertex; they are pair- wise shared by the cubes /, II, V and VI. The correspondence 
between the links viewed from the viewpoint of the cubes //, V, VI with those of the cube 
/ is given by the following table: 
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II 


I 


1 


7 


16 


6 


18 


10 


b 


g 



V 


I 


9 


11 


16 


4 


17 


3 


b 


d 



VI 


I 


2 


8 


17 


5 


18 


12 


b 


e 



Performing the same steps as in deriving the 6j symbol for the vertex a we arrive to the 
following result for the vertex b: 




(112) 



We notice that vertex a is of the 'even' and vertex b is of the 'odd' type: all other vertices 
of the lattice can be considered as either 'even' or 'odd'. Therefore, eqs. ( |111| , |112|) give 
actually the full result. Combining them together we find that a sign factor 



-1) 



-v 



(113) 



should be attributed to all links of the lattice. 

Let us prove that this sign factor is equivalent (in the vacuum!) to a sign factor 



(-1) 2J = (-1) 



-23 



(114) 



attributed to all plaquettes of the lattice. We recall that all links are shared by two even 
cubes whose faces carry plaquette values J. We first attribute all links to only one (out of 
the two possible) cubes, according to some rule. Many such rules can be suggested, the only 
requirement being that each link is attributed to one and only one even cube. An example 
is given by the following construction: we choose the edges 12,5,9,2 and 7 (see Fig. 2) as 
'belonging' to the cube shown on that figure. The rest six edges will then 'belong' to one of 
the neighbouring even cubes. For example, the edge 1 will be counted as 'belonging' to the 
cube II (see Fig.3). Indeed, from the cube II point of view that edge will be of the type 7, 
and so forth. It can be seen that, in these scheme, every link of the full lattice will 'belong' 
to one and only one even cube. 
We have, therefore, a sign factor 



\2ji2+2j 5 +2j9+2j 2 +2.j7 



;ii5) 



attributed to the cube I. Next, we recall that, e.g., jxi enters the 3jm symbol together with 
the plaquette angular momenta Jb and Je (see (|9~C|)). Using eq. ( |9"2"D appropriate to the case 
we can replace (-1) 2 -? 12 = (-l) 2J s+2J B> Similarly, (-1) 2 -? 5 = (_i)2J s +2Jf ; an d so on As a 
result we get that the sign factor ( |115| ) is equal to 

/_^\2J A +2J B +2J C +2J D +2J E +2J F (116) 

This procedure can be repeated for all even cubes of the lattice. It proves the above statement 
that the product of all link sign factors ( |113| ) can be replaced by the product of all plaquette 
sign factors (|114|) . It should be stressed that this proof is valid only for the vacuum, i.e. for 
the partition function itself but, generally speaking, not for the averages of operators. 
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Appendix D. 9j symbols from the Wilson loop 



Let the Wilson loop in the representation j s go through the links ...,15,1,17,..., see Fig. 3 for 
notations. It means that one has now to integrate three D-functions of the link variables 
f 15,1,... j instead of two, as it was in eqs. ( |103| , |108| ) of the previous Appendix, the rest 
integrations remaining unchanged. We have now 



Jl Jl Js 

o a u a m a 



Jl Jl Js 

Ob u b m b 



117) 



^ 15 ^ ) _ 9e (^ 5 )^| ai _ Se (C/ 1 t 5 )^ me (C/ 15 

i\m e -m a ( Jl5 Jl5 Js \ ( Jl5 J15 js 



'1181 



Using the other 3jm symbols related to the vertex a (see eqs.( |99|]T02|) ) and the Kronecker 
symbols from eqs. ( |1 04j - |l 07| ) we get for the vertex a: 



"a" 



^ ^ -J y — s+tu+p— m 



Jl J4 Jl2 

o r z 



Ji J13 Ju 
—u —w —p 



JU JVo Jl2 

p q z 



J4 J15 Jl3 

r s w 



Jl Jl Js 

o u m 



J15 Jib Js 
q s m 



(119) 



[we notice that r + z = —o, w + p = —u and that o + u + m = 0, hence the sign factor is +1; 
we change the signs of all projections in the second 3jm symbol, and permute the columns 
in other 3jm symbols to match the definition of the 9j symbols as given by eq. 



J4 Ji J12 

(_l)ii-i4+il4+il«+* <j jg 

jl3 j[ JU 

To get the final sign factor we have used the relation (— I) ±2 ii ±2 i2±2i3 
.71,2,3 originating from one 3jm symbol. 

Acting in the same fashion we obtain for the vertex b: 



(120) 

+ 1 valid for any 



i \ —v—y—t—x+m 



Jl J9 J2 

—o —y —v 



Jib Jis Ji 
t x u 



Jig Jn J9 
t q y 



Jis J2 Jn 

X V s 



Jl Jl }s 

o u m 



J2 Jl J9 
.^[-h+he+hr+js I fij jg jl7 

jl8 j'l jl6 



Jn Jn Js 
q s m 



;i2i) 
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